Rising bosonic electroweak virtual effects 
at high energy e + e colliders 

M. Beccaria 

Dipartimento di Fisica, Universita di Lecce and INFN, Sezione di Lecce 
Via Arnesano, 73100 Lecce, Italy 



G. Montagna and F. Piccinini 

Dipartimento di Fisica Nucleare e Teorica, Universita di Pavia and INFN, Sezione di Pavia 

Via A. Bassi 6, 27100 Pavia, Italy 



OO 



(N 



- 1—1 

x 



F.M. Renard 

Physique Mathematique et Theorique, UMR 5825 
Universite Montpellier II, F-34095 Montpellier Cedex 5. 



C. Verzegnassi 

Dipartimento di Fisica, Universita di Lecce and INFN, Sezione di Lecce 
Via Arnesano, 73100 Lecce, Italy 
(February 1, 2008) 



The virtual massive gauge boson effects on several observables in the process of charged fermion 
pairs production at e + e~ colliders are computed up to the TeV energy range in the Standard Model, 
working in the "Z-peak subtracted" representation. It is shown that these effects increase regularly 
with energy, approaching the typical ten percent relative size. A careful numerical evaluation is 
proposed, where the conditions dictated by gauge invariance are consistently implemented. 
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OO ■ I. INTRODUCTION 

The accurate computation of virtual electroweak effects has appeared to be a fundamental theoretical contribution 
Oh! f° r & correct interpretation of the experimental results obtained in several years of running at the Z resonance at 
LEP1 and SLC Q. Leaving aside the remarkable degree of accuracy to which the Standard Model (SM) has been 
\ tested at the one loop level, the prediction for the value of the top mass (spectacularly confirmed by its experimental 
discovery) and the subsequent limits on the "low" value of the Higgs mass undoubtedly represent an illustration of 
^ ' how relevant a rigorous calculation of virtual effects may be. 

From a purely technical point of view, a calculation of electroweak (i.e. of not conventional QED type) virtual 
effects appears to be somehow simplified on top of the Z resonance, since for trivial reasons the non resonant terms 
can be evidently neglected. This eliminates the full set of "box" diagrams, thus reducing the theoretical task to that 
of computing self-energies and vertices. In the self-energies, the dominant effect is provided by light fermion pairs 
that determine the running of a^ ED ^ up to the Z mass [|J, and by the heavy top that shifts from 1 the value of the 
p parameter || . 

A welcome feature of these fermionic contributions is that they are separately gauge independent. This is also true 
for what concerns the remaining "large" ~ ml virtual effect provided by the Zbb vertex [Q. Therefore, the full set 
of "relevant" virtual electroweak contributions at the Z peak appears to be relatively "simple" , the main role being 
played by either the fermionic content of self-energies or by a gauge independent subset of the bosonic effects. 

As soon as one moves away from the special LEP1/SLC configuration, the expectation is that the previous simple 
features might disappear. In particular, the role of boxes (where s-channel W exchange is intrinsically not gauge 
independent), could become much more relevant. Should this be the case, one extra immediate complication would 
arise since the gauge dependence e.g. of W bosons must be cancelled, in the contributions to physical observables, 
by that of suitable combinations of corresponding self-energies and vertices where virtual Ws are exchanged. Clearly, 
only the sum of such special box, self-energies and vertices combinations can acquire a physical meaning, and must 
be computed as carefully as possible. 

To make the previous discussion more quantitative, we can now provide one numerical example, i.e. the calculation 
of the box contribution to cross sections and asymmetries in the LEP2 energy range. Working in the t'Hooft £ = 1 
gauge, one sees || that, at the largest considered energy y/s ~ 200 GeV, the related effect can become quite respectable 
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(a relative one percent in the muon cross section and a two-three percent in the hadronic cross section). Since this 
effect is apparently increasing with energy, a question that naturally arises at this stage is that of whether this feature 
will remain valid, or stronger, at the future linear e + e~ colliders, whose energies are expected to lie in the 500 GeV 
range (or more j6|). In particular, the problem that one should face is that of computing correctly the relevant gauge- 
invariant combinations of self-energies, vertices and boxes. In fact, if their size became "quite" respectable, the extra 
problem of the calculation of the initial-state QED radiation corrections should also be considered in a self-consistent 
way i.e. including all the various terms in the convolution process that is usually involved |Q. 

The aim of this paper is precisely that of computing, in a fully gauge invariant way, the contribution of bosonic 
effects to physical observables of e + e _ annihilation into charged fermion pairs. This will be done in an energy range 
from 200 GeV to 1 TeV (although we could increase the upper limit if requested), working systematically in our chosen 
subtracted theoretical framework. As we shall see, the numerical effect generally increases with energy, approaching 
a typical ten percent size at the TeV scale. This will be fully discussed in the following Sec. || and Sec. [II, where the 



numerical calculation of a few relevant observables will be performed without taking photon emission into account. 
The large size of the effect will then motivate a rigorous calculation of the initial-state QED radiation corrections; 
this will be treated in Sec. |v]. A final brief discussion will then be given in the concluding Sec. 0. An appendix 
contains explicit formulae that are involved in the procedure but not essential for a qualitative understanding of the 
main results. 

II. BOSONIC VIRTUAL EFFECTS IN THE "Z-PEAK SUBTRACTED" REPRESENTATION 

It has been recently shown |^,|^] that the process of e + e~ annihilation into fermion pairs, at squared center of 
mass (cm.) energy q 2 larger than M§ , can be conveniently described using an effective parameterization, strictly 
valid at the one loop level, where only subtracted quantities appear as "electroweak corrections" . In particular, the 
contributions to all observables generated by the Z boson exchange are systematically subtracted on top of the Z 
peak, at q 2 = M§ , in full analogy with the photon contributions that are subtracted "on top" of the photon peak at 
q 2 = 0. As a consequence of this procedure, the theoretical input in the various Z contributions will contain quantities 
(partial Z widths and asymmetries) measured at the Z resonance. These input parameters will replace the commonly 
used Fermi coupling G M , and we devote the interested reader to Refs. ]||| for an exhaustive discussion of the details 
of the approach. The final result is that, in the complete expression of a general observable of the process e + e~ — > //, 
there will be some "leading" terms containing ckqedi Mz and Z peak observables and four independent one loop 
subtracted corrections that contain all the residual physical information and are defined as: 

K,e f {q 2 -e) = F 77 , e /(0, 9) - F 77 , e/ (g 2 , 9) (2.1) 
R ef (q 2 ,9) = Iz, ef (q 2 ,9) - Iz,e f {M 2 ,9) (2.2) 

\H Z (n 2 0\ - ^Z,ef{q 2 ,9) A lZ ,ef(M 2 ,9) 

yZ-rf„2 a\ _ A Zl , e f(q 2 ,9) A Zl , e f{Mj,9) 



Iz,ef{q\6) = 2 q — r2 [F zz , ef (q 2 ,9) - F zz>e/ (Mj|, 0)] (2.5) 



In the previous expressions, we have used the following definitions 

q 2 
q 2 - M 

Azz, ef {q 2 ,9)= A z 2^(0,9) +q 2 F Z z,ef{q 2 , 9) (2.6) 

Azz, ef (q 2 ,9) = A zz (q 2 ) - (q 2 M$)[(T$,v}$) + (tJJ,^) + (q 2 - M 2 )A^ f (q 2 ,9)} (2.7) 

F^ ef (q 2 ,9) = F^ ef {q 2 ) - (r« „$) _ (T™,v™) - <?A%$tf,9) (2.8) 
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A 7 z,ef(q 2 ,0) 



A^ z (q 2 ) g 2 -M| 

„2 V „2 



)( r /Z/'°/- / 



(? 2 - M z) A \z°ef(<l 2 ' ) 



(2.9) 



q 2 q 2 q 2 



-Cr$, «$) - (« 2 - ^i)4te1(9 2 ^) (2.io) 

The quantities j4y(g 2 ) = Ay(0) + q 2 Fij(q 2 ) = 7, if) are the conventional transverse 7,Z self-energies. 

^77°7Z Z7 zz e/ (*? 2 ' ^) are ^ ne projections on the photon and Z Lorentz structures of the box contributions to the 
scattering amplitude ^4 e y and the various brackets (r^, u„) are the projections of the vertices on the different Lorentz 
structures to which A 77 , Azz, A 7 z, Az^ belong. In our notations A^ ox ^ is the component of the scattering amplitude 
at one loop that appears in the form v^lA^^v^^ where v^f = — |eo|Q e ,/7 M is what we call the photon Lorentz 
structure and analogous definitions are obtained for A^ x \ , A^° x ^ with the Z Lorentz structure defined 



as v 



( Z )<V - _M_ n ,u./„ „0 „,5 



■j e j = — 2s co ^' J '^Ve f ~ 9AejJ 5 )- We do not insist here on these technical points that are, we repeat, fully 
discussed in previous references • Our notations are almost rigorously following the conventions of Degrassi and 

Sirlin Jl^,[l3|], and our approach is often motivated by t heir p revi ous observations. In particular, one can immediately 
understand that all the four quantities listed in Eqs. ( |2.l| )-( [j~4| ) are separately gauge-independent, since they con- 
tribute different Lorentz structures of the process. As a consequence of this fact, whenever computing contributions 
e.g. from virtual Ws exchanges, that are intrinsically gauge dependent, it will be essential for a clean interpretation 
of the result to stick the various effects in the proper previous combinations. 




(a) (b) (c) 

FIG. 1. Feynman diagrams for self-energy WW components. $ and are the would-be Goldstone bosons and ghosts, 
respectively. 



To illustrate the previous qualitative statements with a concrete example, we start from the calculation of the 
charged boson contributions to the generalized photonic correction A Qje /(q 2 , 9) of Eq. (2.1). Since this quantity is by 



construction gauge-independent, we have performed the calculation in the familiar t'Hooft gauge. This means that 
the contributions from the unphysical charged would-be Goldstone bosons must be incorporated and added to those 
coming from the charged Ws (and ghosts). We shall follow the order given by the various definitions and compute 
the transverse self-energies terms first. In the case of A Qe j, these correspond to the Feynman diagrams shown in Fig. 
1. A straightforward calculation leads to the result (only the sum of all separate terms of Fig. 1 is given): 

K%! ww <(i 2 ) = f£ +1,+1 '(o) - f^'+'V) = 
£){jf*H-£<l II -|r l(1 -' T)l 

I OU2 pi 2 

^ ' 11 ' dx In li-J^l-s)!} (2.11) 
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FIG. 2. Feynman diagrams for WW vertex component (Fig. 2a), for single W vertex (Fig. 2b) and for fermion self-energy 
diagrams (Fig. 2c). 



The previous self-energy component A*' 



(WW) 



ef 



is universal i.e. independent of the final fermion features. This 



is not true for the remaining vertex and box contributions. It is convenient therefore to consider the various cases 
separately by treating final muons (or taus), final u quarks and final d quarks states in succession. Starting with final 
muons, we have first to consider the diagram involving Ws in Fig. 2a (we are assuming massless final fermions, so that 
their couplings with would-be Goldstone <& bosons vanish). There is also no photon-neutrino coupling so no diagram 
of the type of Fig. 2b. In principle one should finally add the fermion self-energy diagrams of Fig. 2c. However 
their contribution is ^-independent and disappears in the subtraction. This is also the case for any type of tadpole 
diagram, that we therefore never discuss here. This leads to the following expression of the (two Ws) contribution to 

the generalized vertex r M 7 (g 2 ), to be conventionally defined as T^ W \q 2 ): 



r (WW)(2a+2b) ^ _ 

+2 / dxidx 2 [3 In \M w (xi + x 2 ) - q 2 X\X 2 



7 M (l- 7 5 ){-2(--7) + 



lnM w 



q 2 (xi + x 2 - xix 2 ) 
(M w (xi + x 2 ) - q 2 X\X 2 ) 



(2.12) 



where e = (4 — n)/2 and 7 = ln(47r) + 7, 7 being the Eulero-Mascheroni constant. 

Projecting this vertex on the photon Lorentz structure, according to our procedure, gives then 



AS W (9 2 ) = -(-) / d Xl dx 2 [5 



q 2 X\X 2 



q 2 (x 1 



x 2 - XiX 2 ) 



M 2 v {x l 



x 2 ) 



M w { Xl 



+ x 2 ) - q 2 xix 2 



(2.13) 



For a final uu or dd pair, we have to include in the calculation of r„ ~(q ) another (single W) diagram as shown 
in Fig. 2b. This is, in fact, the non universal component of such charged boson vertex, since to the graph of Fig. 2a 
with two Ws a universal quantity is associated. Strictly speaking, diagrams with charged would-be Goldstone bosons 
should also be included when the mass of the involved quarks cannot be neglected. This corresponds to the cases of 
final ti pairs and of final bb pairs. Top production will not be considered in this paper since our technique, based 
on Z peak subtraction, does no t ap ply in this case; for bb final states the specific calculation of that contribution to 
the vertex will be given in Sec. |A l| . For "light" u, d quarks we then obtain the following contributions (we give the 
projection on the A aie f(q 2 ,0) combination directly): 



A vert.(W)( 2 



a 
Gtt 



dx\dx 2 [In |1 



q x\x 2 



.(1 



X1-X2) 



q 2 (l- Xl )(l- X2 ) 
r(l — xx — x 2 ) — q 2 x\x 2 



(2.14) 



A 



vert.(W) 



(q 2 ) = 2AZ^Hq 2 



(2.15) 





FIG. 3. Feynman diagrams for WW box contribution; Isf = — 5 (Fig. 3a), -^3/ = +5 (Fig. 3b). 
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The final step is now the addition of the box contributions, that correspond to the single Feynman diagram of Fig. 3 
((a) for final states with I^f — — i, (b) for I 3 f = +i) in the cases where the mass of the final quarks can be neglected 
(in practice, for our purposes, the only exception would be the box with a final bb pair, to be discussed in the appendix 
together with the related vertex) . The expression of the WW box contribution to the full scattering amplitude has 
been already given in the literature ]l4j . By a straightforward projection on the photon Lorentz structures one easily 
obtains for final ^t + /i~, uu, dd pairs 

Af « {WW) {q\9) = - V e ) 2 [I(q 2 , t, M W ) + I 5 (q 2 ,t, M w )} (2.16) 

Af % {WW \q\ 9) = J^L(1 - « e )(l - v u )[I(q 2 ,u, M W ) - I 5 (q 2 ,u, M w )} (2.17) 

KZ (WW \q\9) = 3^4(1 - «e)(l - « d )[/(g 2 ,t,M w ) + I 5 (q 2 ,t,M w )] (2.18) 
where = 1 — 4|Q/|s 2 and I(q 2 ,t, Mw), h{q 2 , t, My/),I[q 2 ,u, Mw), h{q 2 , u , Mw) are known functions [] defined in 



Sec. A3, with t = -%(1 - cos6), u = + cos 



The full charged boson contribution is then given by the sum 

^charged W ){q ^ e) = ^(WW) {q2) + (WW ){q2) + £«ert. (W) (g2) + £B« (V U , [q2 f) } ( , [t)) 



in which the self-energy component is universal, the vertex with two Ws is also universal except for the m 2 dependent 
terms for bb final states and the vertex with one W and the box terms depend on the final state. 

The calculation of the charged boson contribution to the three remaining gauge independent combinations R, V lZ , 
V Zl proceeds now in strict analogy with that illustrated for A Q e y, replacing the initial photon with an initial Z . The 
calculations are straightforward and do not r equi re a spec ial mention. Rather than illustrating them in detail with a 
long series of formulae, we just refer to Sees. |A l| and |A 3| in which all the individual terms of the following sums are 



given: 



R (charged W) ^ = ( WW) (?2 j + R vert. ( WW) ((?2) + R vert. (W) {q2) 

-.Be 



R?. ox(WW) (q 2 ,6) (2.20) 



yTZ (charged W) ^ = y'f (WW) ^ + y^Z, vert. (WW) ^ + y-fZ, vert. (W) 

+Vjf' Box{WW \q 2 ,e) (2.21) 



yZj (charged W) ^ _ ys.e. (WW) / 2\ + yZ-y, vert. (WW) ,2^ + yZ-y, vert.(W) / 2\ 



+V^- B "' ww \q 2 ,e) (2.22) 



From the previous expressions Eqs. (pT9|)-(p32j), we are now in a position to derive the charged boson contribution 



to any observable of the process e + e~ — > // (f ^ e). This is made possible by the expression that the related 
differential cross section assumes in our formalism. We shall be limited for the moment to the situation where 
polarized electron beams are not considered (this might be a very interesting case in future measurements, as stressed 
in a previous paper |l5f| ). In this case, we shall write: 

%q 2 , 9) = |(1 + cos 2 Q)a\ S + cosOaf (2.23) 



dcosO ' 8 

where 



1 We thank W. Hollik for sending us his Fortran implementation of these functions. 
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[l + 2A Q , e/ ((7 2 ,l 



-i-WrV\in w q2 ~ M % i[ J1 ni/2[ -' U/2 

1 V((<7 2 - Mf ) 2 + M|r|) U M z J L JV/(Mf)Mz (1 + v^(l + g»)V 



r 3r eil 



3r, 



:[l + A Qie/ ( g 2 ! 0))-i? e/ (g 2 



3T, 



+ 



Wlz il N f (M^)Mz i 

( q 2 - m 2 ) 2 + m 2 t 2 z 



(2.24) 



^ = ^(¥){M0)IQ/|[ 



( g 2 -Af 2 ) u 3r 1/2 3iy 

<?W - M 2 )2 + M 2 T 2 z ) l[ M z l L JV>(Ml)M ; 



ll/2 



[ ][ 

' ^^[l + A a . ef (q 2 ,6)~R ef ( q 2 ,6)}+ - w " 



3r e ir 3Tf 



(1 + 5 2 ) 1 /2(1 + ^2 ) 1/ 

(l+v 2 )(l+^.) v e ej 



( q 2 -M 2 ) 2 + M 2 T 2 z 

m v e f( q \e)}}} 



(2.25) 



where Nf(q 2 ) is the c onve nt ional color factor which contains standard QCD corrections at variable <j 2 and the theo- 
retical input of Eqs. (2. 24), (2. 25) contains T e , Tf, v e = 1 — 4s 2 , Vf = 1 — 4|Q/|s 2 measured on top of Z resonance 
with an accuracy that is sufficient to prevent any dangerous theoretical uncertainty in our predictions for future e + e~ 
colliders, given their expect ed re a listic experimental accuracy, as exhaustively discussed in Refs. ^[|9|. 

From the previous Eqs. J2~19|)-(|2~2l) it is now perfectly clear how to compute, in a rigorously gauge-independent 
way, the charged W contribution to various observables. This is the result of the insertion in those equations of the 
charged boson contributions that we have listed. 




H H 
FIG. 4. Feynman diagrams for neutral boson contributions to Z self-energy. H and $ are the Higgs boson and the would-be 

Goldstone boson, respectively. 




G 




In a quite analogous way, one would compute the overall neutral boson component or the fermionic contribution to 
only self-energies, that is separately gauge independent. In fact, we have also evaluated the neutral boson contribution 
by an analogous calculation, involving Z self-energies (Fig. 4), Z vertices (Fig. 5) and ZZ boxes (Fig. 6) (including 
when necessary the neutral would-be Goldstone bosons and Higgs exchanges). Without entering the details of the 
procedure, we anticipate that these terms contribute, for reasonable values of the Higgs mass, all observables by 
an amount that is much smaller than that coming from the charged boson diagrams. The light fermion self-energy 
contribution (Fig. 7) is, on the contrary, relevant and must be suitably taken into account for a complete numerical 
estimate. This will be discussed in detail in the two forthcoming Sections. 



/ 




FIG. 7. Feynman diagram for fermionic contribution to self-energy correction. 



III. PREDICTIONS FOR OBSERVABLES IN ABSENCE OF PHOTON EMISSION 



Our final step has been at this point that of translating all the previously illustrated calculations into realistic 
numerical predictions for the various observables of the process e + e~ — > //. We have concentrated our attention 
in an energy range \fq* between the Z resonance and 1 TeV, although we could have considered higher energies as 
well. Technically speaking, we have divided our programs into two separate phases. In the first one, we have written 
analytic expressions for the four quantities A Q , R, V lZ , V Zl that should contain the bulk of the numerically relevant 
components (we have systematically ignored, as one might have noticed, imaginary parts in the various quantities; 
this approximation should be quite tolerable at high energies). More precisely, we have taken into account in a 
conventional way, i.e. including leading-log resummation, the running of QqeDi due to fermion loops, that appears 
in A Q (and plays a fundamental role) Q and we have also inserted the full fermion contribution to the Z and 7Z 
self-energies. These terms, added to the leading terms that contain a, Mz and the Z-peak quantities, will constitute 
what we shall call "non bosonic component" of the various observables. To these expressions we have added those 
given in Eqs. (2.19)-( |2.22 ), plus the negligible neutral boson component. These terms will generate the charged and 
neutral boson contributions to all observables. 

The second phase of our approach has been that of implementing our theoretical expressions into a semi-analytical 
program called PALM (PAviaLecceMontpellier program). This program, which is already available, will be fully 
illustrated in a separate dedicated publication. In PALM all the integrations over the Feynman variables yielding the 
one loop subtracted self-energy and vertex corrections are performed numerically [fl6| according to semi-analytical 
algorithms optimized to guarantee reliable and accurate results for the electroweak effects. In particular, the two- 
dimensional integrals providing the vertex contributions are first analytically reduced to one-dimensional integrals, 



2 The leptonic contribution to oqed is exactly included as well as the top-quark effect, both below and above the tt threshold. 
The light quark contribution has been implemented following the generally adopted parameterization of Eidelman-Jegerlehner 
in Ref. [01. 
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that are then treated numerically on the same ground as the self-energy corrections. A numerical integration over 
the fermion scattering angle is performed for weak boxes in order to get their contribution to the integrated cross 
section and forward-backward asymmetry. It has been checked that the numerical procedure adopted provides results 
in perfect agreement with the formulae obtained by analytic integration over the Feynman parameters. The complete 
one loop electroweak contributions to A a ,ef, Ref, V^f and V^ 1 are codified in PALM following a completely modular 
structure, where the corrections contributing each Lorentz structure are grouped into a single routine. This will allow 
a quite easy implementation in the program of analogous virtual effects, such those due to SUSY and other new 
physics models. 

To check the validity of our computation we have first compared our numerical results, without taking QED 
radiation into account, with the corresponding ones of the most recent version of the program TOPAZO 11711. The 
input parameters used for the comparison are Mz = 91.1867 GeV, m t0 p = 160 GeV, rrim ggs = 115 GeV and 
a s (Mz) = 0.120, together with the most recent data for Z peak observables needed as input in PALM fljsfl . For the 
sake of comparison, the value mt op = 160 GeV has been chosen because fits to all electroweak data except the direct 
determination of m t0 p and A% prefer a "low" top mass and therefore this choice turns out to be consistent with 
the input Z peak observables of PALM. Further, for a meaningful comparison with the TOPAZO predictions, the input 
Z partial widths of PALM have been corrected for the effect of standard QED final-state correction. For simplicity 
we have not included at this level in our program the aforementioned heavy top contributions to the b vertex and 
box. These actually only contribute the bb differential cross section, and will affect the total hadronic quantities that 
we have considered (in practice, at this preliminary stage, the cross section for production of the five light u, d, s, c, 
b quarks) to a presumably negligible extent (we shall return to this point later on). 



Relative dev. for ct 5 
Relative dev. for rj, 
Absolute dev. for A FF 




V s (GeV) 

FIG. 8. Comparison between the predictions of PALM and those of TOPAZO for the leptonic and hadronic cross section oi 
(I 7^ e) and 05 (relative deviations) and for the leptonic forward-backward asymmetry Aps,i (absolute deviation). 



The comparison has been carried through in an energy range yq^ between Mz and 1 TeV. It is depicted graphically 
in Fig. 8 for three observables, chosen on the basis of their expected good experimental accuracy: the cross section 07 
for production of charged leptons (I ^ e), the related forward-backward asymmetry Afb,i and the cross section for 
five light quark production 05. 

As one sees from Fig. 8, the agreement between our calculations and those of TOPAZO is indeed impressive in 
the full considered energy range. This is particularly true for the two leptonic observables, for which the difference 
between the two approaches is always less than about two permille, well below the expected experimental accuracy, 
estimated to be of about a relative one percent |l9| . This is also true for the hadronic cross section, with the exception 
of two small regions around 300 GeV and at the end of the range when approaches 1 TeV, where the relative 
difference reaches a five per mille value that is, though, still below the one percent expected experimental accuracy. 
In fact, we believe that this discrepancy can be originated by the heavy top effect in the b vertex and box that we did 
not include in our program, and are contained in TOPAZO. Since this difference is indeed modest, we feel that it can 
be safely neglected at this stage. Note that we did not consider the intrinsic theoretical errors implicitly contained 
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in both programs (e.g. PALM has as theoretical input a number of Z widths and asymmetries). We estimated such 
uncertainties to be, typically, of the permille size in both cases, which would still reduce the (already miserable) 
existing discrepancy. 

Having checked in this way the reliability of our calculation, we now move to the real goal of this paper, that is the 
consideration of the (charged) boson effects. These are shown for the three separate variables tr^, Apg,^ an d 05 in 
Fig. 9 as computed by our program PALM. 

As one clearly sees from Fig. 9 the effect of such bosonic contributions becomes indeed relevant as soon as the cm. 
energy raises above a typical value of approximately 200 GeV, that corresponds to the LEP2 limit. At energies of 
about 500 GeV, in the expected range of a future e + e~ linear collider (LC) the relative size of the bosonic terms 
is already larger than two percent in and six percent in 175. This relative effect continues to rise when the cm. 
energy increases, reaching values of the ten percent size in the hadronic cross section in the 1 TeV limit. This clearly 
generates, in our opinion, two independent and rather relevant theoretical issues. 



t> ; Without ISR 

0.08 - 
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. To:a gauge invariant effect 
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FIG. 9. The rising charged boson effects as functions of the cm. energy, for the three observables cr M , 05 and Afb,^, in 
absence of photon emission. The separate effect of charged boson self-energies and vertices (dash-dotted line) , WW box (dashed 
line) and their gauge invariant sum (solid line) is plotted. The relative effect of charged boson effects is shown with respect to 
the observables including the gauge independent overall neutral boson component and fermionic corrections. 

The first problem is that of understanding the reasons of such a visible rise. At first sight, this appears to be 
stronger than a purely logarithmic one, as it can be easily checked by numerical inspection. Should this be indeed the 
case, one could not advocate the expected asymptotic purely logarithmic increase dictated by renormalization group 
arguments that would apply to the transverse self-energies (in their generalized gauge independent definition proposed 







in Refs. [ pL 2|Jl 3| , that also includes a small amount of "pinched" vertex ]20[| ). The extra increase should therefore be 
rather addressed to the remaining component of the various combinations (by this we mean what would remain after 
the purely logarithmic, separately gauge independent, "generalized" Degrassi-Sirlin self-energies have been isolated). 

As a matter of fact, it has been pointed out in a previous reference [Ell, and stressed again in Ref. |2fl, that the 
non pinch part of the W vertex, Fig. 2b, behaves in the large q 2 limit as the squared logarithm of q^M^- Such 
a behaviour is often denoted as a "Sudakov logarithm" [^2| since the limit q 2 /M^ — > oo, is formally identical to 
an "infrared", QED like, Mw — ► limit. On general QED analogy arguments, one would expect such terms to 
appear also in the box contributions to the various combinations. More practically, one actually sees such terms in 
the analytic expressions of boxes given in Sec. A3. Thus, the existence of such a behaviour at large q 2 appears to us 
theoretically expected. 

Strictly related with this increase is, in our opinion, the problem of how to treat perturbatively such effects when 
they become too large. In other words, one should study a mechanism of rcsummation analogous to that proposed 
by Sudakov |2^] in the QED case. Although this problem is at the moment of rather academic interest (since after 
all even at 1 TeV the effect is still below ten percent), it might become rather relevant in the case of future colliders 
at energies in the few TeV range. 

The second problem that arises is a consequence of the fact that the bosonic effect is, in any case, sizeable when one 
moves beyond the LEP2 energy limit. This means that a proper estimate of photon emission effects, that also includes 
such full bosonic terms in the quantities to be convoluted, becomes imperative. We shall perform this calculation in 
the next Sec. IV . This will require a technical discussion which we tried to make as concise as possible. 
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FIG. 10. The same as Fig. 9, including QED initial-state radiation (ISR). The cut \fp > 100 GeV, s' being the invariant 
mass of the event after ISR, is imposed to exclude the Z radiative return. 
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IV. TREATMENT OF PHOTON EMISSION EFFECTS 



For the aim of the present study, we are interested in the contribution of the initial-state radiation (ISR), since it 
is well known that it constitutes the largely dominating QED effect.^] The evaluation of QED radiation effects has 
been obtained in our calculation adopting the QED Structure Function (SF) method [Q. ft provides a simple, non 
perturbative approach, based on factorization theorems of infrared and collinear singularities, for computing leading- 
log (LL) photonic corrections to any arbitrary e + e~ hard-scattering cross section. The method allows to keep under 
control the long-distance universal logarithmic contributions, due to the soft and/or collinear photon radiation, by 
convoluting a kernel cross section with electron(positron) structure functions. For a typical event selection where a 
cut on the invariant mass of the event after ISR is imposed, the QED corrected cross section can be simply calculated 



as 



a(s) 



1 — Xcu 



dx H(x, s)<7q ((1 — x)s) 



(4.1) 



where H (x, s) is the radiator or flux function, and <r is the elementary de-convoluted cross section eventually including, 
as in our case, the process-dependent pure weak and QCD corrections. The radiator represents the probability that 
a fraction x of the whole cm. energy s is carried away by ISR and is calculable as: 



H{x,s) 



dz 



D(z,s)D 



l-x 



(4.2) 



where D(x,s) is the electron SF. A typical solution for the radiator, generally used for LEP1/SLC calculations, 
reads ^,p3|: 

H(x, s) = A2PX - 1 + hi(x, s) + h 2 (x, s) 
hi(x,s) = --(3(2- x) 



h 2 (x,s) = -I3 2 



lnfl — x) 

(2 - x) (3 ln(l - x) - 41nx) - 4— - + x - 6 



/3 = 2(^)[L-1], L = ln( S /m 2 ), A 2 = 1 + + (^)' 



5 2 . 



(4.3) 



The first exponentiated term in H(x,s) describes all order soft multi-photon emission, hx^XjS) are associated to 
hard bremsstrahlung in collinear approximation, (3 being the QED c ollin ear factor. The perturbative coefficient A2 is 
the 0(a 2 ) soft+virtual K factor for s-channel processes. Equation (4.3) for the QED radiator, that accurately works 
for precision studies of the electroweak interactions on top of the Z resonance, becomes inadequate far from it, where, 
as discussed in Ref. |p4| , higher-order hard photon effects may become numerically relevant as a consequence of the 
phenomenon of Z radiative return. Actually as pointed out in Ref. p4| , 0(a 2 ) next-to-leading (NL) hard-photon 
corrections, that are known from the exact 0(a 2 ) diagrammatic calculation ]2f| and numerically negligible close to 
the Z peak, introduce measurable effects at LEP2 and beyond and can be taken into account with the replacement 



H(x, s) — > H(x, s) — h,2(x, s) + $2 (1 — x, s) 



(4.4) 



where the explicit expression for 82 can be found in Ref. j2f|. These 0(a 2 L) corrections are of the same order of 
magnitude as the third-order LL ones, that come from the emission of three hard collinear photons and have been 
recently analyzed in the literature for their impact on two-fermion production cross sections at LEP2 |24| . Including 
both NL 0(a 2 ) and LL O(0 3 ) contributions, an accurate QED radiator in analytic form can be cast as follows VM: 



3 In principle, one should also consider other conventional QED effects as those due to final-state radiation and initial-final 
state interference. However, for the experimental set up considered in the following, the (very small) effect of final-state 
radiation is automatically included in our approach by means of the input physical quantities V f. Further, the contribution of 
initial-final state interference and QED boxes has been evaluated by using the previously quoted program TOPAZO and found 
to be negligible as compared to the effect of ISR. 
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H(x, s) = A 3 f3x p - 1 + hi(x, s) + <5f (1 -x,s) + h 3 {x, s) 



2T 1 5 1 

— + — a; + 4(1- -x) (tt 2 - 6 In 2 x + 3Li 2 (x)) 
6 3 \ , 2 , J _ . 4 . 7 



+31n(l -x) 7 -x +ln (1 - x) -7+ - + -x 

x 2 / V x 2 



A 3 = A 2 



-6 lnx(6 

r a \ 3 
7T 



x) + 6 In x ln(l — x) ( 6 3x 



(L-lf ( A-i7T 2 --^ 2 )(l) 

v ' 1 16 2 3^ w 



(4.5) 



where i[)( n >(z) is the ri-th order p olyg amma function, ip^ n \z) — d n ip(z)/dz n , ip( z ) = r'(z)/r(z). Equation (4.1), 
together with the radiator of Eq. (4.5), is the formulation implemented in PALM for the simulation of initial-state 
photonic effects to the muon and hadronic cross section. 

Since the numerator of the forward-backward asymmetry is a less inclusive quantity than a fully integrated cross 
section, the convolution integral for the asymmetry is different from that of the total cross section and is given by |27[| 



A FB {s) = 



1 



a(s) 



4(1 — x) 

dx - H FB {x, s) a FB ({l - x)s) 
(2 x) 



H FB {x, s) = H{x, s) + h%£Ax, s) + h FB 2 (x, s) 



(4.6) 



where a FB (s) is the forward minus backward cross section and a(s) is the total cross section as given by Eq. (|4.l|). 
The radiator H FB (x, s) can be written as 



(4.7) 



hpg x stands for 0(a) NL corrections to s-channel processes |25 



uLL 

IL FBa 



where H (x, s) is the radiator of Eq. ([i 
0(a 2 ) LL "antisymmetric" corrections [^9|. Equation (4.6), with the radiator of Eq. (4.7), is the convolution integral 
implemented in PALM for the calculation of QE D ra diatio n eff ects to the muon asymmetry. 

In the program PALM the integrations of Eqs. (4.1) and (4.6) are performed numerically by using a semi-analytical 
procedure improved by standard importance sampling tricks to take care of the infrared and Z radiative return 
peaking behaviour. 

The impact of the initial-state photon emission on rising electroweak virtual effects is shown in Fig. 10 and Fig. 11, 
for two different typically adopted selection criteria, i.e. excluding the Z return (y/P > 100 GeV, with s' = (1 — x)s) 
and including it (y/s* > 20 GeV). As it can be clearly seen, a proper inclusion of QED radiation is unavoidable to 
get a reliable evaluation of bosonic virtual effects. Actually, when excluding the Z return, the QED convolution can 
give a relative reduction of the order of 25-30% of the overall bosonic effects to and 0-5 and a significant lowering 
of the same effects to A FBill . When including the Z return, the impact of the QED convolution is much larger, 
reaching a 60-80% relative effect in the cross sections and changing sign to the net effect to the forward-backward 
asymmetry. The event selection with the inclusion of Z return, in particular, is responsible of a much flatter behaviour 
of the bosonic corrections, whereas the exclusion of Z radiative return reduces the absolute size of the corrections 
without critically modifying their peculiar rising slope. Typically, the size of the effect on the hadronic cross section 
at 500 GeV is reduced by QED convolution from six percent to four percent, which certainly represents a drastic 
difference visible at an expected experimental accuracy below the 1% level. 
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FIG. 11. The same as Fig. 10, with VV > 20 GeV, to include the Z radiative return. 



V. CONCLUDING REMARKS 



We have shown in this paper that an accurate computation of theoretical predictions for the process of e + e~ 
annihilation into charged fermion pairs at future high energy colliders requires the full determination of the bosonic 
contributions, including boxes which are usually neglected on top of the Z resonance. We have performed the relevant 
calculations in the t'Hooft gauge, but the overall contributions that we have shown are completely gauge-independent. 
For this overall effect we conclude that, since the contribution from its "generalized" gauge independent self-energies 
can only be of logarithmic type when q 2 becomes large, extra stronger effects of Sudakov type can only be produced 
by the combination of the "non pinched" parts of vertices with the boxes. These rising effects are evident in our 
calculation, and become visible as soon as one crosses the 200 GeV range. As we have shown, the effect of a proper 
inclusion of QED radiation is indeed numerically essential and should be carefully taken into account in the theoretical 
predictions. 

One possible interesting generalization of our investigation would be the search for analogous sizeable contributions 
in a SUSY calculation of virtual effects at future colliders. On a qualitative basis, we would expect the presence of 
important effects produced by either chargino or even neutralino vertices and boxes pof . Such a possibility is at the 
moment under investigation. 
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APPENDIX: COMPLETE 1-LOOP ELECTROWEAK CONTRIBUTIONS TO A a , EF , Ref, AND 



We give here the explicit expressions of the quantities that appear in Eqs. (2.19)-(2.22) 



1. Charged boson contributions 

Universal self-energy WW components (Fig. 1): 



+ ^ + ^ ) r &in i i -^ (:i -" 11 



,„ „ /' t 2x(l-x)[20c 4 + (l-2 S 2 ) 2 }-20c\ q 2 

He f [q) ~ ( ^>J X{ 4^ [ q 2 -Ml L 

x{l~x)M 2 2 -17 + 36s 2 -16s 4 g 2 - 2 

M 2 , - x(l x)M 2 J + 4c 2 , 2 \ 2 - Mf £ ~ Mz£zJ} 



- In 1 J^Z " ^ I 



-rln 1 - -Vxfl - ar) =■ In 1 

n 2 M 2 M 2 M 2 

q lvl w Z 1V1 W 



KJ- <"'"'> (rt = -(f) /"<fa4l— »4' - 5 V + ±i^£r(l - x),'] In |1 - J£-*(1 - *)| 

JO Q cs S CS 

. 1 Nr 8s 2 -9 2 r c 2 11 -12s 2 ,„ | M 2 ,„ 

M§ cs " 5 ; M| + ^s— " X)Mz] ln 11 ^ Mf X{1 ~ 

WW vertex component (Fig. 2a, universal for light fermions) 

~ AV er,(WW ){q2) ^_ { a } f f dxidx2[3lnll _ 9^2 I i ? 2 (-l+-2-^2) , 



e/ W; ir y s 2 >JJ x l M^(x 1 +x 2 )-M^x 1 x 2 l 

| (g 2 - M|)M 2 v ,(xi + x 2 )(xi + .t 2 - xix 2 ) 
[M 2 v {x 1 + X2 )-q 2 x 1 x 2 ][M 2 v {x 1 +x 2 )-M 2 x 1 x 2 ] 



} 
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(<? 2 - M|)M^(xi + x 2 ){xi +x 2 - xix 2 ) 



+ 



[M w ( Xl + x 2 ) - q 2 x 1 x 2 ][M 2 v (x 1 + x 2 ) - M\x x x 2 \ 



q 2 - M 2 Z q 2 x lX2 , q 2 (x 1 +x 2 -x 1 x 2 ) 



H^H^ vertex component for final state 

T^ert. (WW),„2x _« f f, , „ , ™t M 2 , ( Xl + X 2 ) + TTl 2 (1 - gi - X 2 ) 



~2~^ 



(A9) 



<f' """'V) = £5 / /dT,dr 2 {[12c 2 + (1 - 2» 2 ) 

M w (x! + z 2 ) + m 2 (l - xi - x 2 ) - g 2 xix 2 i . 2r 9 2 ( x i + x 2 - ZlZ 2 ) 



_m|_ 

M 2 J 



M^,(ii + x 2 ) + m 2 (l - xi - x 2 ) - M 2 z x x x 2 M^(x\ + x 2 ) + m 2 (l - x\ - x 2 ) - q 2 x x x 2 

M z (x\ + x 2 - x\x 2 ) 



M w (xi + x 2 ) + m 2 (l - xi - x 2 ) - M z xix 2 2 



}+ l R vert.(WW ){q2) (An) 



~fZ,vert. (WW), 2\ _ a ( C \ f f , , r,, | M w {x X + X 2 ) - q 2 X\X 2 



Vr Ve "- V ""V) = £0 I ld Xl dx 2 {3\n\ 
+ 



w >-2^s>J J -^^^M 2 w {x 1+ x 2 )-M 2 zXl x^ 
(q 2 - M 2 )M 2 v {x 1 + x 2 )( Xl +x 2 - Xl x 2 ) 



[M w ( Xl + x 2 ) - q 2 x 1 x 2 ][M w (x 1 + x 2 ) - M 2 x lX2 ] 
q 2 -M z ac ff M 2 v {x 1 +x 2 )+m 2 {l-x l -x 2 )-q 2 x l x 2 
+ ^JJ rfXlrfX2[3lD ' M 2 w{ x 1+ x 2 ) +m 2 { l- Xl -x 2 ) 1 

+77T7 ^+ 2 n~ XlX2) > 2 ]> (A12) 

M 2 / (x 1 +x 2 )+w 2 {\-x 1 -x 2 )-q 2 x 1 x 2 s 



v Z^er, (WW ){q2)= a_ { ^ Jj dxidx2 



,q 2 -M 2 q 2 x lX2 q 2 (xi + x 2 - x x x 2 ) 



q 2 M^( Xl +x 2 ) M$y(x! +x 2 ) - q 2 Xl x 2 

2 



+ ( 8^ } / / dx l dx 2 {[12c 2 + {l-2s 2 ) 

. M w (xi + x 2 ) + m 2 (l - xi - x 2 ) - q 2 xix 2 



ml 1 



' M w( x i + x ?) + m l(^ -xx- x 2 ) - M\x\x 2 

+4c 2r q 2 (xi +X 2 - XlX 2 ) M z (xx +X 2 - XjX 2 ) 

M w( x i + x ?) + m t( l - x i - x 2) - q 2 x\x 2 M w (xi + x 2 ) + m 2 t (l - xi - x 2 ) - M 2 z x x x 2 

(A13) 

Single W vertex component (Fig. 2b, light fcrmions; non universal) 



KZt iW \q\0) = ^ (' [ d Xl dx 2 [In |1 - . f XlX2 -| - n , 2 f, {1 ~ Xl)i \~ X2 l ] (A14) 

6tt J J M w (l- Xl -x 2 y M 2 v {l-x l -x 2 )-q 2 x 1 x 2 i 
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2irs 2 J J M^(\ - x x - x 2 ) - q 2 x x x 2 

, (g 2 - MpM^(l - Xl - x 2 ){l - Xl ){\ - x 2 ) 



[M^(l - Xl - x 2 ) - g 2 iii 2 ][M2,(l - an - x 2 ) - Mf x l3 ; 2 ] 



(M/) (? 2 ) = (i-y)i?r- (H/ V) 



kt {w \<i 2 ) = ^-\)Kr- (H V) 



} 



KT' <"V) = ^[[ d^dx 2{ ln\ M J^- X ^^\~ M t lX2 \ 
M 4ttcs J J M^(l - xi - x 2 ) - q 2 x lX2 

| (g 2 -M|)M 2 / (l-x 1 -x 2 )(l- a;i )(l- a;2 ) 

[M^(l - xi - x 2 ) - g 2 xix 2 ][M2 / (l - xi - x 2 ) - M§xix 2 ] 

v2 Z, V er, (W ){q 2 j9) = yver, {W) ^ J J 

g 2 xix 2 | g 2 (l - xi)(l - x 2 ) 

ikf2,(l - xi - x 2 ) M2 / (l-Xi-x 2 )-q 2 XiX 2 l 

v: z, ver , {W){q2 ^ = VeT , (w){q2) (t^Ml^ 1 1 dxidx2{ 

r, ,-. q 2 x lX2 q 2 (l - xi)(l - x 2 ) 

1111 M2,(l-xi-x 2 )' M 2 v {\-x 1 -x 2 )-q 2 x 1 x 2 ) 

v Zy,vert. m { ^ 0) = (1 _ lf!) K -*- W(g 2 ) 
y Z 7lW ert. W(g 2 )(9) = (1 _ ^f!) Kr *. (W) {q 2 } 

Single W vertex component for bb final state 

Tvert. (W), 2s _ « , m t x /* /* . . r,_ I ^jU 1 ~ X l ~ gg) + W 2 (x X + X 2 ) - g 2 XiX 2 

A - i6 (9 , -3^ ll+ 2M^ , yy 1 2{ 1 M 2 / (l-xi-x 2 )+m 2 (xi+x 2 ) 
m 2 m t 2 + q 2 (l - xi)(l - x 2 ) 



M^(l - xi - x 2 ) + m 2 (xi + x 2 ) M^r(l - xi - x 2 ) + m 2 (xi + x 2 ) - g 2 xix 2 ' 

4s 2 N 4s 2 m 2 

2~ 

W 



jyvert. (W), 2x « I A A Hon ^ \ ^ " l t 

R * {q)= 8^ dXldX2mi -— ] -^W 



[-ln 
+ 



Myy(l - xi - x 2 ) + m 2 (xi + x 2 ) - q 2 x\x 2 
M 2 / (l - xi - x 2 ) + to 2 (xi + x 2 ) - M|xix 2 

g 2 (l-xi)(l-x 2 ) M|(l-xi)(l-x 2 ) 



M 2 (1 - xi - x 2 ) + m 2 (xi + x 2 ) - g 2 xix 2 M 2 (1 - X\ - x 2 ) + m 2 (xi + x 2 ) - Mfxix 2 



3 3 J tl M^(l-xi-x 2 )+m 2 (xi+x 2 )-g 2 xix 2 



M^(l -xi -x 2 ) +m 2 (xi +x 2 ) - M|xiX 2 JJ 2 e ^ 
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, n Z,vert. (WO, 2x a ffj j n \ M w( l ~ X ± ~ X 2) ~ M Z x l x 2 \ 

Aires J J ' M 2 v (l-x 1 -x 2 )-q 2 x 1 x 2 ' 

(g 2 - M 2 )M^(1 - x x - x 2 )(l - gi)(l - gg) i 
[M^(l - xi - x 2 ) - g 2 xia; 2 ][M^(l - xi - x 2 ) - MfxiXa] 1 
ac q 2 - M 2 Z m\ /*/*,, r , , - Xi - x 2 ) + m 2 (x 1 + Xg) - g 2 XiX 2 

+ ^ } y y dxidx2h ln 1 - X1 - , 2) + ^ + X2 ) 1 

+ m 2 + q 2 (l-x 1 )(l-x 2 ) m| 

M^(l -xi- x 2 ) +m 2 (xi +x 2 ) - q 2 xix 2 M^(l - xi - x 2 ) + m 2 (xi + x 2 ) 

T ,Z7,uert. (W), 2x , « uro/1 4s 2 8s 2 ™ 2 /• /■ 

(<z) = ( w ){[2(1 -— } u^ ] JJ dxidX2 

M$y(l - Xi - x 2 ) + m 2 (xi + x 2 ) - g 2 xix 2 
n M^(l - xi - x 2 ) + m 2 (xi + x 2 ) - M|xix 2 

+ g 2 (l-x 1 )(l-x 2 ) M 2 (l- Xl )(l-x 2 ) 

M^(l - xi - x 2 ) + m 2 (xi + x 2 ) - g 2 xix 2 M^(l - x\ - x 2 ) + m 2 {x\ + x 2 ) - Mf xix 2 

_ 4s 2 _mf_ _ 8s^ f f 2 1 

+ K 3 3 7y dXldX2mt[ M* r {l-x 1 -x 2 )+m?{x 1 +x 2 )-q*x 1 x 2 

\ 1 (A27) 

M 2 ,(l - xi - x 2 ) + m 2 (x! + x 2 ) - M 2 X!X 2 J 



2. Neutral boson effects 

ZH self-energy (Fig. 4) 



^• e ' (Z "V) = (g^) jf rfx {[M 2 (x - 2) + M|(l - *)][^^|« 

x(l-x) 1 Mfx(l-x) n 

+ n/r2„2 1 iT72 T\ 77 + 775" ln I 1 — 



M|x 2 + M 2 r (l-x) M| 1 Mfx + M 2 (1 - x) IJ 

M 2 .x + M^(l-.x)- g 2 x(l-x) 
W = ln| M 2 x 2 +M^(l-x) 1 (A29) 

W-iuii q 2 x(l-x) 1 Mfx(l - x) 



Single Z triangle (Fig. 5) 



- ^V) = ^ ,^ " dx ld x 2 [-ln\l 



yH >~ 32ns 2 c 2 J J |X Af 2 (1 - X! - x 2 ) 1 

g 2 (l-x 1 )(l-x 2 ) 

+ M|(l-x 1 -x 2 )- g 2 x 1 x 2 J 1 J 

ef ' 32?rs 2 c 2 J J L 1 Mf(l - xi - x 2 - xix 2 ) 1 
q 2 {l-x 1 ){l-x 2 ) (1 -x^l -x 2 ) 

M§(1 — X! — X 2 ) — 9 2 X!X 2 1 — Xi — X 2 — XiX 2 
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„ 7 Z,vert. (Z) , 2 , ,0^eQ- ~ vj) ^ f [, , r , , M|(l ~ ^1 ~ X 2 ) - g^i^z 
V ef ) = ( C ,„„3„3 ) / / ^ldx 2 [-ln|— 



e ^ v 647rs 3 c 3 ^y y 1 'i i _ Sl _ a 2 _ arixa) 

g 2 (l-a;i)(l-x 2 ) (l-s^l-aa) ^ 

M|(l — a;i — ^2) — g 2 a;ia;2 1 — xi — 22 — a:i:r2 
f-MZa\Q f \v f ) f f dxidx2[ _ lnll t^l 



q 2 Attsc J J M\(\ — x\—X2) 

q 2 (l- Xl ){l-x 2 ) , 



M|(l ~ xi — X2) — q 2 x\x 2 ' 



y M\Qf\ns*c3' J J 21 1 M*(l-x 1 -x 2 -x 1 x 2 ) 

q 2 (l-x 1 )(l-x 2 ) (l- Xl )(l-x 2 ) 



Mg(l — x\ — x 2 ) — q 2 x\x 2 1 — X\ — x 2 — x\X 2 
q 2 -M 2 z _^ ) f [ dxidx2[ _ ln]1 _ fax. 



q 2 A-ksc J J M\(\ — x\—x 2 ) 

q 2 {\-x 1 ){\-x 2 ) 



M§(1 — X\ — x 2 ) — q 2 X\X 2 



(A33) 



(A34) 



3. Box diagrams 

In the light fermionic case the WW and ZZ box amplitudes have been written in Ref. juj in terms of I and J5 
functions. We have also written them in the Passarino-Veltman description |3l|| in terms of Di functions as defined in 
Ref. |3^] . We have checked |33| that both expressions agree numerically using the package of Ref. j34j] . The extension 
to the 66 case could be done very easily by just putting m = wit instead of m = inside Di(0, 0, 0, 0, t, q 2 ; M, 0, m, M). 

Projecting these amplitudes on photon and Z structures then give the following results. 

WW box (Fig. 3) 

Defining from Fig. 3a 

n w n w_< n t [I(q 2 ,t,M w ) + I 5 (q 2 ,t,M w )} 
u = D d = [D 27 - -(D 25 - D 2i - D u )\ = — (A35) 



and from Fig. 3b (D having the meaning of changing t into it inside D) 



D™ = [AD 27 + \{Diz + 2D 26) + 1 -{D 13 + 2D 25) + |(D 1X + D 12 + 2 Am)] 

= [I(q 2 ,u,M w )-I 5 (q 2 ,u,M w )} 
2q 2 



we obtain: 



(A36) 



X*? WW) iM = l^ku^--e)(l-v f )DY (A37) 



R i.Bo X ,WW ){q2e) = c 2 {M 2 -q 2 )a D y 



v - t z (B OX ,WW ){ 2 e) = c(M 2 - g 2 )(l - v e )a 



yZl (Bo^ww) 2 = c(M z -q )(l-v f )a 
3 47r|<5/|s ci J 
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ZZ box (Fig. 6) 

For the ZZ box, we obtain 



with 



A£^ ZZ V > 0) = -( 25 J Q g ' g4c4 )[(l - - e 2 )(l - t^)Df + v e v f (l - vl){l v))Dl] (A41) 

Ri Bo X ,ZZ ){q ^ e) = _ h <Mz-jf) [AVeV{D Z + + vl)[l + v})D Z 2] (A42) 

v ,z (B OX ,zz ){q ^ e) = -h f a ^~p \2v f {l - v*)D? + v e (l vl)(l + v))D z 2 ] (A43) 

v g (B OX ,zz ){q ^ e) = -( a Wz-/J )[2v e (l vj)Df + v f (l + vl)(l v).)D z 2 ] (A44) 

Df = [D 27 + ^^-(D 24 - D 25 + D u ] 

-[4^27 + y (Dia + 2£> 2 6) + -(Dia + 2£> 25) + - (.Du + D12 + 2D 24 )] 



[/(g 2 , i,M z ) -/(<?>, M z ) 



2g 2 



(A45) 



Dl = (<^ 25 - D 24 - D U) ] = [W,^) + 7 5 ( g 2 ,,,M,)] (A46) 



4. Fermion pair component 



The contribution from Fig. 7 is universal and leads to (to/ being the mass of the fermion in the loop, gvf = v fhf, 
9 as =hf)- 

£«.e. (//)( g 2 ) = 2^/ / 1 dx2a;(l-.x)ln|l-^a;(l-a;)| (A47) 

7T m f 



z(l-a;) 1 . M_ 



77T ln |l--^(l-^ 



mj - M%x(l - x) M| to/ 

V-Af 2 m 2 -M§a;(l-a;) J 



-2(<£/ + - (A48) 



to? — g 2 ai(l — a;) 

1 o 2 1 M 2 

J = -2 ln I 1 - T^C 1 - *)! - 772 ln I 1 - tJF *(i - *)l ( A5 °) 



<7 2 m 2 M| to 2 . 
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